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ABSTRACT 

Aims. The main goal is to study the dynamics of the gravitationally stratified, field-free cavities in the solar atmosphere, located under 
small-scale, cylindrical magnetic canopies, in response to explosive events in the lower-lying regions (due to granulation, small-scale 
magnetic reconnection, etc.). 

Methods. We derive the two-dimensional Klein-Gordon equation for isothermal density perturbations in cylindrical coordinates. The 
equation is first solved by a standard normal mode analysis in order to obtain the free oscillation spectrum of the cavity. Then, the 
equation is solved in the case of impulsive forcing associated to a pressure pulse specified in the lower-lying regions. 
Results. The normal mode analysis shows that the entire cylindrical cavity of granular dimensions tends to oscillate with frequencies 
of 5-8 mHz and also with the atmospheric cut-off frequency. Furthermore, the passage of a pressure pulse, excited in the convection 
zone, sets up a wake in the cavity oscillating with the same cut-off frequency. The wake oscillations can resonate with the free 
oscillation modes, which leads to an enhanced observed oscillation power. 

Conclusions. The resonant oscillations of these cavities explain the observed power halos near magnetic network cores and active 
regions. 
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1. Introduction 

Observations show a high spectral power of oscillations (v > 
5 mHz) in the vicinity of active regions. The observed velocity 
' power maps are sometimes referred to as 'photospheric power 
halos' (Braun et al. 1992; Brown et al. 1992; Hindman & Brown 
1998; Jain & Haber 2002; Muglach et al. 2005; Moretti et al. 
2007; Nagashima et al. 2007; Hanasoge 2008). The power en- 
hancement is observed only in the Doppler velocity power maps 
and not in continuum intensity (Hindman & Brown 1998; Jain & 
Haber 2002; Muglach et al. 2005; Nagashima et al. 2007). The 
enhanced power spectra peak at the period of ~ 3 min. 

In addition, observations show that the chromospheric ac- 
tive regions and magnetic network elements are suiTounded by 
"magnetic shadows", which lack oscillatory power in the higher 
frequency range (Mcintosh and Judge 2001; Krijger et al. 2001; 
Vecchio et al. 2007). From the observations it seems reasonable 
to conclude that both the photospheric power halos and the chro- 
mospheric magnetic shadows reflect the same physical process. 

In our previous paper (Kuridze et al. 2008), we proposed an 
explanation of the photospheric power halos as the effect of a 
magnetic canopy on the wave dynamics. We showed that the 
field-free cavity regions under the magnetic canopy can trap 
high-frequency acoustic oscillations, leading to the observed 
increased high-frequency power in the photosphere, while the 
lower-frequency oscillations are channeled upwards in the form 
of magneto-acoustic waves (Erdelyi et al. 2007; Srivastava et al. 
2008). However, those calculations have been performed with- 



out taking into account the gravitational stratification, which is 
important at the photospheric level. 

In the present paper, we study acoustic oscillations in the 
stratified, cylindrical field-free cavity regions under the mag- 
netic canopy as a possible explanation of the power halo phe- 
nomenon. There are various analytical and numerical investi- 
gations of the acoustic wave propagation in isothermal strati- 
fied atmospheres (e.g.. Lamb 1908; Rae & Roberts 1982; Fleck 
& Schmitz 1991; Kalkofen et al. 1993; Sutmann et al. 1998; 
Roberts 2004). However, almost all these calculations involve 
simple ID models (for vertically propagating acoustic waves), 
because 2D and 3D models necessarily include internal grav- 
ity waves, which considerably complicates the analysis. In or- 
der to avoid further complications due to the gravity waves, it 
is possible, however, to only consider isothermal propagation, 
which automatically neglects gravity waves. This enables us to 
study the 2D propagation of only isothermal acoustic waves or 
pulses in a stratified atmosphere. Here, we study normal isother- 
mal acoustic oscillations and the propagation of acoustic pulses 
in a cylindrical field-free cavity by solving the 2D Klein-Gordon 
equation in cylindrical geometry. Unlike the ID case, the two- 
dimensional solution allows us to investigate the variation of the 
density/velocity perturbation amplitude in the azimuthal direc- 
tion. First, we use a standard normal mode analysis to study the 
spectrum of possible acoustic oscillations in the cavity region. 
Then, we solve the Klein-Gordon equation with impulsive forc- 
ing in the photospheric region. 
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The outline of the paper is as follows. In Sect. 2 we present 
the basic hydrodynamic equations for the problem of the propa- 
gation of acoustic waves in a stratified cavity medium. Section 3 
describes the normal mode approach and the resulting oscilla- 
tion spectrum in the cavity. In Sections 5 and 6 we analyze the 
response of the cavity region to a pressure pulse. The results ob- 
tained are discussed in Sect. 7. 



2. The wave equation 

We use the ideal hydrodynamic equations for a gravitationaUy 
stratified field-free cavity: 



|.V.(pv) = 0, 
p(|+v.Vv) 



-Vp+pg, 



dp „ „ „ 
+y - Vp + ypV • V = 0, 

ot 



(1) 
(2) 
(3) 



where v is the fluid velocity, p and p denote the pressure and 
density, respectively, y is the ratio of specific heats, and g is the 
gravitational acceleration. 

The considered acoustic cavity under the cylindrical mag- 
netic canopy has a semicircular form (Fig.l). Therefore, it is 
convenient to consider the equations (1-3) in a cylindrical co- 
ordinate system. For simplicity, we study the 2D case, i.e. per- 
turbations polarized in the (r, (p) plane, where r is the radial co- 
ordinate and <p is the azimuthal angle. 

In this cyUndrical (or polar) coordinate system, the equilib- 
rium pressure and density are related as foUows: 



dp( 



'01 



dr 

1 dpoi 
r d<l> 



= grPOl, 



= g^POl, 



(4) 



(5) 



where poi> Poi denote flie equilibrium pressure and density and 

gr = -gsincp, 
= -gcos4>, 

are flie r- and (^-components of the gravitational acceleration 
vector, g, respectively. 

The Eqs. (4-5) then easily lead to the equilibrium values of 
the pressure and the density in the cavity: 



Poi - Poie ^ , 

« % simp/A 

Poi =Poie ^' , 
where 
A=-^, 

is the pressure scale hight. The atmosphere is considered to be 
isothermal. Therefore, this pressure scale height is constant. 

In the considered cyUndrical (polar) coordinate system, the 
Unearized hydrodynamic equations for adiabatic fluctuations can 
be written as: 



(6) 
(7) 

(8) 
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dt 
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dpi 
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r dip 



■g cos (fipu 



dpi odpi Cq 2 

— Cg— + Poi —(sin (pUr + cos ipu^)a)., = 0, 

ot ot g 



(10) 
(11) 
(12) 



where Cg = ypoi/poi is adiabatic uniform sound speed and a)b 
denotes the Brunt- Viiisala frequency: 



(01 = 



( y-i)g 

yA 



(13) 



These equations describe the acoustic and internal grav- 
ity waves propagating in a stratified atmosphere. In principle, 
file cavity under a magnetic canopy may trap both types of 
waves. Unfortunately, the simultaneous consideration of both 
wave types is very compUcated, especially from an analytical 
point of view. However, by adopting specific approximations, 
we can consider each of the wave types separately. The limit 
of incompressibility, e.g., neglects the acoustic branch and only 
gravity waves remain in this case. On the other hand, the limit of 
isothermal propagation neglects the gravity branch of the spec- 
trum and in this case only the acoustic branch remains. The lat- 
ter limit means that the temperature exchange is so rapid that 
any temperature fluctuation in the perturbations is zero. This is 
achieved when y approaches unity. 

We use the approximation/hmit of isothermal propagation in 
the remaining part of this paper. For isothermal perturbations, 
Eqs. (9-12) lead to the following wave equation: 



,5^1 
'°dr^ 



502 



rsin^c + — — -H 

r ) or r'- 



, & dpi 

+ - COS(p—— 

r dip 



d'pi 

dt^ 



(14) 



Equation (14) can be written in a more convenient form by 
using the following transformation: 



pi =pie 



-/Irsin^ 



(15) 



where the parameter A can be chosen in such a way that the 
resulting wave equation has the form of a Klein-Gordon equa- 
tion in the cylindrical/polar coordinate system. This procedure 
is analogous to the weU-known transformation in a ID Cartesian 
coordinate system (Kahn 1990; Sutmann et al. 1998). 
By substituting Eq. (15) into Eq. (14) we obtain: 



2^^Pl 



^0 dpi ^ cl d^pi 



d% 

dr^ " r dr " dtp^ dt^ 



{2Acl 



\ ■ dpi 
■g)smip— 



-{2Acl-gy-^^-^^(A'cl-Ag)pi(r)^0, 

from which we find the foUowing condition on A 
2A4-g = 0, 



or 



dr 



dip 



2c^ 



1 

2A' 



(16) 



(17) 



(18) 
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Fig. 1. Simple schematic picture of small-scale magnetic canopy 
overlying a field-free cavity with gravitational stratification, g. 
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Using this condition, we obtain 



r dr 



1 d% 



1 



where (J^^ refers to the acoustic cut-off frequency 



4A 



4A2' 



(19) 



(20) 



Equation (19) is the Klein-Gordon equation in cylindri- 
cal/polar coordinates. This equation governs the isothermal 
propagation of acoustic waves (and pulses) in the considered 
gravitationally stratified medium. The very dynamic convective 
layer located under the cavity can excite quasi-harmonic wave 
trains as well as acoustic pulses. Clearly, the cavity will re- 
spond in different ways to these different types of perturbations. 
We study two different responses of the cavity to such different 
perturbations from below. First, we consider the cavity as res- 
onator for acoustic waves and, using the standard normal mode 
approach, we study the spectrum of possible wave harmonics 
trapped in the canopy. Next, we consider the propagation of a 
pressure pulse, which can be caused in the lower regions e.g. 
due to the eruption of new granular cells or by some small-scale 
magnetic reconnection processes. 

3. Normal mode analysis 

A standard Fourier analysis of Eq. (19) with respect to time, 
leads to the equation: 



j^^pi dpi d'^px 2 



+ r h 



+ r 



,2 \ 



Pi 



0. 



which gives 



We consider the solution of Eq. (21) in three different cases 
depending on the size of the circular frequency w with respect to 



When CO - ojac, equation (21) leads to the Laplace equation 
in cylindrical coordinates. The solution of this equation inside 
the canopy (i.e. inside the semicircle r - rc, where is the 



Fig. 2. The oscillation frequency, w, versus the cavity radius, rc, 
for the first zeros of the m = 1 (bottom line), m-2 (middle line) 
and m-Ti (top line) harmonics. 



cavity/canopy boundary (see Fig. 1)) with the condition p\ - 
pc((f>) (i.e. a Dirichlet condition) along the boundary is given by 
(Morse & Feshbach 1953) 



pi - ^[A„, cos(m0) -i- B„, sin(m0)] 



where 

Am - 

and 



2 r 

Pc(a, 

^ Jo 



) cos(ma)da. 



B, 



2 r. , 

Pcia 

^ Jo 



) sin(ma)Ja. 



(22) 



(23) 



(24) 



This means that one of the free oscillation modes in the cavity 
oscillates at the cut-off frequency for any m (which plays the role 
of azimuthal wave number). For the photospheric sound speed 
of 7.5 km/s, the acoustic cut-off frequency is of the order of 
0.031 s^^ with a corresponding period of ~ 200 s. 

In the case o) > lOac, Eq. (21), again after a Fourier analysis 
with respect to the azimuthal direction and with azimuthal wave 
number m, transforms to a Bessel equation with the general so- 
lution 



(21) ^ ciJ„(kir) + C2Y„(kir), 



(25) 



where 7„, and Y„, are the bessel functions of the first and second 
kind, respectively, and ci and C2 are arbitrary constants. 

When cj < ojac, then a similar Fourier analysis transforms 
Eq. (21) to the modified Bessel equation, with the general solu- 
tion 



pi = cn,I,n(kir) + C4K„(kir), 



(26) 



where /,„ and K,,, are the modified Bessel functions of the first 
and second kind, respectively, and cj and q are arbitrary con- 
stants. The density perturbations must be finite at r = 0, which 
yields the condition C2 - C4 = 0. We also need a second condi- 
tion at the upper boundary of the cavity in order to find a unique 
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solution of equation (21). The correct approach is to find an an- 
alytical solution of the magnetohydrodynamic equations in the 
overlying magnetic canopy region and then merge it with the so- 
lution (25) at the canopy/cavity interface (r = r^). This would 
yield the exact oscillation spectrum of the cavity, which may 
reveal some wave leakage. Unfortunately, we could not get an 
analytical solution of the complete cylindrical MHD equations 
for the gravitationally stratified plasma in the canopy region. 
However, we can exploit the fact that the density quickly de- 
creases from the photosphere to the chromosphere. According 
to the VAL-C atmospheric model (Vernazza et al. 1981), there is 
the sharp density gradient in the low area of the chromosphere 
(Fontenla et al. 1990). Therefore, it is justified to approximate 
the sharp gradient between the photospheric and the chromo- 
spheric densities as a discontinuity at the canopy/cavity inter- 
face. The density is then much higher in the cavity than in the 
overlying magnetic canopy (Fig. 1). This allows the use of a free 
boundary condition, i.e. pi = 0, at the cavity/canopy interface, 
which yields an approximate spectrum of acoustic oscillations in 
the cavity. The influence of overlying canopy can be studied by 
complete numerical simulations, but this is not the scope of the 
present paper and should be done in the future. The free bound- 
ary condition leads to the expressions: 

ciJmihrc) = 0, for oj > (Oac, (27) 

c^Imihrc) = 0, for oj < ojac, (28) 

The first condition gives oscillatory solutions (oscillating with 
respect to the r- coordinate), while the second condition is sat- 
isfied only when C3 = 0. Therefore, Eq. (28) yields the trivial 
solution of Eq. (21), and in the following we will concentrate on 
the oscillatory solution connected with the condition (27) only. 

The spectrum of acoustic oscillations in the cavity can then 
be deduced from the zeros of the solution Jm{z), which (as it is 
known from the theory of Bessel functions) are all real when m > 
-1, and they can be easily found in tables (see e.g. Abramowitz 
& Stegun 1967). 

Figure 2 shows the dependence of the frequencies of the 
acoustic oscillations on the cavity radius, r^, for the first zero 
of Jmiz) for the w - 1 , 2, 3 harmonics. It is evident from this 
figure that the frequency decreases with increasing r^, as can 
be expected from physical considerations. This means that these 
cavities, with typical granular radii = 400 - 800 km, support 
oscillations with frequencies in the range 0.07 - 0.05 s"' (i.e. 
with periods of 1.5 - 2 min). For considerably large cavity radii, 
the frequencies of all the harmonics tend to the cut-ofl" value. 

It seems that cylindrical cavities support higher frequency 
oscillations than typical quiet Sun regions. This is because 
in cylindrical cavities, the oscillations are partly along the 
^-direction due to the small ki term for the first zero of Bessel 
function in Eq. (27), which prevents their leakage upwards and 
leads to their trapping in the cavities. On the other hand, the 
regions without overlying canopies (e.g., nonmagnetic quiet 
Sun regions outside the network cores) can not trap the high 
frequency oscillations, because they may propagate upwards. 
Therefore, it is expected that the high frequency oscillations are 
trapped around the magnetic network cores and active regions 
where the cylindrical magnetic canopies can be formed. The pe- 
riods of the lower order harmonics in cavities with a radius of 
1000 km (Fig. 2), correspond to the observed spectrum of acous- 
tic oscillations (5-8 mHz) (Moretti et al. 2007) (see flieir Fig. 1). 

The free boundary condition, which is used to obtain the os- 
cillation spectrum, does not lead to wave leakage in the canopy 




400 600 800 1000 1200 400 500 600 700 800 900 1000 
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Fig. 3. Real cOr {m = \ bottom Une, m = 2 top line on the left 

panel) and imaginary w, (m = 1 top line, m - 2 bottom line on 
the right panel) parts of wave frequency vs the radius of field- 
free cavity r^. The straight line on the left panel is the solution 
with the cut-off frequency. 

region. However, we may explore the problem using the wave 
propagation along a narrow sector around the vertical {<p ^ nl2) 
direction. 

Let us first consider the unperturbed current-free cyUndrical 
magnetic field in the canopy region expressed as (see Kuridze 
et al. 2006) - B^oirdr), where B^q is the magnetic field 
strength at canopy/cavity interface. For simplicity, we consider 
the zero-y8 approximation (Ji = ^npoi/B^ ^ 0), where poi is the 
plasma pressure in the canopy. For this configuration, the solu- 
tion of magnetohydrodynamic equations for the radial velocity 
is found in form of the Hankel function of the first kind (Kuridze 
et al. 2008): 

Ur2 = ic^rH,ni2{k2r), (29) 

where C3 is an arbitrary constant and k2 = w/(2va). Here, 
va = Bj,(}/ ^J4np()2 is the Alfven speed and po2 is the plasma 
mass density in the canopy area. The radial velocity component 
in the narrow sector around the vertical direction of the cavity 
can be obtained combining Eqs. (10), (15) and (25) as 

ftri = e"-'^'^ iclrjkir) + l-J„ihr)) . (30) 

The continuity of the velocity and the total pressure perturba- 
tions at the canopy/cavity interface then leads to the following 
dispersion equation: 

^r/A^^ POl Jmihrc) ^ ^'mjl^^^^c) ^^^^ 

v\ PQ2 J'm^hrc) + Jm(k\ rc)/2A H,n/2(k2rc) ' 

The dispersion relation (31) is a transcendental equation for 
the complex frequency co. A non- vanishing imaginary part of o) 
indicates wave leakage from the field-free cavity into the over- 
lying magnetic canopy. The question then arises how important 
this wave leakage is. The analytical solution of Eq. (31) is com- 
plicated. Therefore, we apply numerical techniques to solve it. 

The numerical solution of Eq. (31) shows that the dispersion 
equation has a real solution with a cut-off frequency o) - oiac, 
and complex solutions with very small imaginary parts. Figure 3 
shows the dependence of the real (left panel) and the imaginary 
(right panel) parts of the frequency on the cavity size r^. The real 
part of the frequency decreases with increasing r^, as expected. 
The damping time of the oscillations = l/w, ~ 10^ - 10^ To 
(ro = 2nloJr) is very large, which indicates there is almost no 
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wave leakage in the canopy region. This means that the first few 
harmonics are trapped in the cavity. It must be mentioned that 
the frequencies of first few harmonics (left panel of the Fig. 3) 
are almost identical to the frequencies obtained with the free- 
boundary condition (Fig. 2). Thus, the oscillation spectrum does 
not depend significantly on the choice of the boundary condi- 
tions. 

It is interesting to find out what happens if the cylindrical 
canopy is replaced by a horizontal magnetic field (this could be 
the case in the quiet Sun regions far from the network cores). 
We have analyzed the spectrum of the field free photospheric 
area under a horizontal magnetic canopy. In this configuration 
there are no trapped high frequency (w > coac) oscillations along 
the vertical direction. The damping times of the vertically prop- 
agating waves are about ~ 2 - 1 To, which indicates there 
leaky nature. On the other hand, horizontally propagating waves 
(with large horizontal wave numbers) have real frequencies. But 
since their velocities are almost horizontal they can not be seen 
in Doppler velocity power maps. Any initial pulse or wave train 
will quickly be dispersed in the horizontal direction. Therefore, 
these modes can not form the observed high-frequency halos. 

4. Excitation by a pressure pulse 

The solar photosphere is very dynamic and it contains many dif- 
ferent types of impulsive sources. For example, the eruption of 
new granules, magnetic field reconnection events just under the 
solar surface, and various other explosive events may take place 
there. Therefore, the excitation of pressure and velocity pulses 
seems to be quiet common under photospheric conditions. In 
this section, we study the propagation of such a pressure pulse 
(which corresponds to a density pulse in the isothermal approx- 
imation) in the cavity region. For simplicity, we assume that the 
pulse has a 5-function shape in space and time. 

Mathematically this corresponds either to the solution of 
Eq. (19) with impulsive initial and boundary conditions or to 
the solution of that equation with an additional term modeling 
the impulsive external forcing. The equation can then be written 
as (with the driving force term in the right-hand side) 



Cavity interface 



Qy.2 f Qy, 



1 d'p. 



1 



"5 Pi ^ 2 ~ 



dfi 



An6{t - to)6(r - ro)6{4> - (po) 



(32) 



where ro, ipo denotes position of the source (Fig. 4). 
Equation (32) governs the propagation of the pulse, which is set 
at r = ro, (p - (po at the time f = to- 
If we consider the case when 



pi(r,t,(p\ro,to,^o) 



dpi{r,(/>, t\ro,to,(po) 
dr 







(33) 



for t < to, then p\ is the Green's function for Eq. (32) and has a 
form (Morse & Feshbach 1953) 



pi(r,f,0|ro,f(),0o) = 



Ao6{t - to - Rico) 



Co 



V(f-fo)2-(i?/co)2 



R 



Hit- to- Rico), (34) 



where An is a constant. 



R 



r-ro 



Irro cos(4) - (po), 




Origin 

Fig. 4. Source point, observation point at f > fo and 
canopy/cavity interface. 



!N 1.5 




Fig. 5. The density fluctuation after the propagation of pulse as 
a function of time and radial coordinate at ^ = 7r/4. The time 
t is normalized by cut-off frequency and r is normalized by the 
photospheric scale height. 

and H denotes the Heaviside step function. 

The first term in expression (34) describes the propagation 
of the pulse. It is followed by a wake which is given by the sec- 
ond term in Eq. (34). This is a demonstration of a well-known 
result (Rae & Roberts 1982; Fleck & Schmitz 1991; Kalkofen 
et al. 1993; Sutmann et al. 1998; Roberts 2004; Zaqarashvili & 
Skhirtladze 2008), viz. the presence of a wavefront which moves 
away from the point (ro, 0o) with the speed cq. The disturbance 
ahead of the wavefront is at rest, but behind the wavefront the 
medium begins to oscillate at the cut-off frequency cogc- 

Using the expressions (15), (18) and (34) the density pertur- 
bation for the wake oscillations takes the form (here fo = is 
chosen without loss of generality): 



Pi 



AoOJacC 



-r sin ip/2A 



- (Rlco)^ 



Co 



Ji - (Rlco)^ Hit - Rico). (35) 



An asymptotic form of the density perturbation can be ob- 
tained using the expression of the Bessel function Ji that holds 
for large arguments it » RIco) (Abramowitz <& Stegun 1967). 
We then obtain 



/TAo-v^e"''™'^^^'^ 
Pi = - V -_ — — cosioJact - 3nl4). 



Co 



,3/2 



(36) 
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After the passage of the pressure pulse, the medium density thus 
oscillates at the cut-off frequency and this oscillation decays 
asymptotically as l/P^^. In fact, the pressure pulse will excite 
the wide spectrum of oscillations, but the higher frequencies har- 
monics escape upwards and the lower frequency harmonics are 
evanescent thus only the oscillations at the cut-off frequency re- 
main. 

It is interesting to study what happens with the velocity 
perturbations in this asymptotic case (t » R/cq). Combining 
Eqs. (10), (11), and (34) we obtain 



/2cos(a.^^.-3;r/4)^,,„^,,^^^ 



dUr AqCq - 

~dt ~ 2Apo 
and 

du^ _ AqCo ^/aJ^ [icosiojact - 3jt/4) 



(37) 



,r sin (p/2A 



The integration of Eqs. (37) and (38) with respect to time and 
using the conditions Ur - and - Q for / — > oo, we obtain 
(see the Appendix for the detailed calculations): 



and 



2 AoVE3^e'-^i°<^/2^cos( 

n pqA fi/2 



cos(a>ac^ - 3;r/4), 



■ COS(L0act - 371/4), 



(39) 



where Ur = Ur/co and = u^/cq. Therefore, the propagation of 
the pressure pulse also sets up a velocity wake oscillating with 
the cut-off frequency. Note, that this velocity oscillation decays 
with time as 1/f i.e. much more slowly than the decay of the 
density wake. Thus the oscillations may easily persist during a 
few tens of cut-off periods (i.e. during ~ 30 - 60 min), which is 
much longer than the expected time interval between consecu- 
tive pulses (which can be estimated as the granular life-time, i.e. 
~ 10 min). Therefore, these wakes of consecutive pulses may 
gradually strengthen the oscillations and thus may lead to the 
observed power halos. 

Figure 5 shows a plot of the density fluctuation, given by 
Eq. (35) (the time is normalized to the cut-off frequency, t' - 
tiOac, and the radial coordinate is normalized by the scale height, 
r' = r/2A). It is seen that the propagation of the pulse is followed 
by a weak amplitude wake which osciUates with the cut-off pe- 
riod. The amplitude of both the pulse and the wake decrease 
with height/time. However, it must be mentioned that the nor- 
malized perturbation of the density, pi/po, increases with r. The 
normalized amplitude of the density perturbation in the asymp- 
totic form, can be written as 



r s'mtp/2A. 



Pi = - 

where 




poco 



Pi 



cos(a)act - 3n/4), 

Pi 
Po' 



(41) 



The normalized amplitude of the asymptotic density pertur- 
bation depends on the direction of propagation (i.e. on the 
0-coordinate). The wake has its maximum amplitude in the ver- 
tical direction (at <p - n/2) and its amplitude decreases for in- 
clined directions (see Fig. 6). 

Figure 7 shows the dependence of the asymptotic veloc- 
ity wake on the azimuthal angle, 0. It is clear that the radial 
component has a maximum amplitude in the vertical direction, 
while the azimuthal velocity component has a maximal ampli- 
tude close to 45". 
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Fig. 6. Normalized asymptotic amplitude of the density wake 
versus the azimuthal angle ^ at r' - 2, /' - 100 (the time t is 
normahzed to the cut-ofF period and r is normalized to the scale 
height). 



(40) 5. The energy propagation 



Our calculations show clearly that att = the energy is concen- 
trated in the initial pressure pulse, and that it propagates away as 
time elapses. The initial pulse, which travels at the sound speed, 
carries most of the energy injected into the medium. We can 
write down an expression for the asymptotic total energy den- 
sity for the wake oscillation, which is the sum of the kinetic and 
acoustic potential energy densities 

sin ip/ A 

E ~ u^ + uI+pI^ kI coi^ {ojact - 37r/4) 



+ Ki 



rs'intp/A. 



(42) 



where. 



and 




PoCo 



It is clear from the expression (42) that the potential part of the 
wake energy decays much faster than the kinetic part. Therefore, 
after some time only the velocity oscillations remain, while the 
density oscillations have decayed. Furthermore, the total energy 
is the largest in the vertical direction (at <p - n/2) and the small- 
est in the horizontal direction (at = 0). It means that the energy 
of the wake oscillations is concentrated near the vertical direc- 
tion (see Fig. 8). Therefore, the wake, which remains after the 
passage of the pressure pulse, tends to oscillate along the verti- 
cal direction. The same result was found by Bodo et al. (1999). 
The total energy of the wake decays in time as ~ l/t. 

6. Discussion and conclusions 

We have carried out a linear perturbation analysis of the two- 
dimensional wave propagation in a stratified, field-free cavity. 
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Fig. 7. Normalized asymptotic amplitudes of the velocity pertur- 
bation components versus at r' = 2, t' = 100 (the time t is 
normalized to the cut-off period and r is normalized to the scale 
height). 



which is located under a chromospheric cyhndrical magnetic 
canopy. We consider the higher density in the cavity than in the 
overlying magnetic canopy region (due to the sharp density gra- 
dient between the photosphere and the chromosphere). Then the 
free boundary condition at canopy/cavity interface leads to the 
spectrum of acoustic oscillation in the cavity. We also have ex- 
amined the response of the cavity to the propagation of a pres- 
sure pulse, which can be excited in the lower regions due to var- 
ious sources, including e.g. small-scale magnetic reconnection 
events. It has been shown that the cavity, which has a radius com- 
parable to the granular dimensions (~ 500 km), supports trapped 
oscillations with periods of 1.5 - 2 min (cf. Fig. 2, 3). However, 
one of the oscillation modes in the cavity is the oscillation with 
the cut-off frequency. Furthermore, the frequency of the acous- 
tic modes, also has the tendency to approach the cut-off value 
for large scale granular size (cf. Fig. 2). On the other hand, it has 
been shown that a pressure pulse propagating through the cav- 
ity sets up a wake, oscillating at the cut-off frequency. The wake 
oscillation will resonate with a free oscillating mode of the cav- 



Fig. 8. Normalized total energy density with respect to <p. It is 
clearly seen that the energy is concentrated near the vertical di- 
rection. 



ity, which may lead to an enhanced acoustic power at the same 
frequency. 

There have been several different explanations of power ha- 
los proposed in the literature: (i) the enhancement of acoustic 
emission by some unknown source (Braun et al. 1992; Brown 
et al. 1992; Jain & Haber 2002), (ii) incompressible oscilla- 
tions, such as Alfven waves or transverse kink waves, in mag- 
netic tubes (Hindman & Brown 1998), (iii) the interaction of 
acoustic waves with the overlying magnetic canopy (Muglach et 
al. 2005; Kuridze et al. 2008) and (iv) a change of the spatial- 
temporal spectrum of the turbulent convection in the magnetic 
field (Jacoutot et al. 2008). However, observations show a lack 
of power halos in the intensity maps (Hindman & Brown 1998; 
Jain & Haber 2002; Muglach et al. 2005; Nagashima et al. 2007), 
which needs an adequate explanation. These observations chal- 
lenge the model of the power halos as due to acoustic waves. 
However, our model easily explains the discrepancy. 

We suggest that the surroundings of the magnetic network 
cores and the active regions consist of many small-scale closed 
magnetic canopy structures (Mcintosh & Judge 2001; Schrijver 
& Title 2003). Granular motions transport the magnetic field 
at the boundaries and consequently create field-free cyhndrical 
cavity areas under the canopy (Fig. 1). The pressure pulses ex- 
cited in lower region propagate through the cavity and leave be- 
hind wakes oscillating at the cut-ofT frequency. These wake os- 
cillations resonate with the free oscillation modes of the cavity. 
Therefore, the cavities accumulate acoustic oscillations with the 
observed frequency. However, the amplitude of the density per- 
turbation in the wakes decreases in time faster (viz. as r^^^) than 
the amplitude of velocity perturbations (viz. as f '). Therefore, 
the density perturbations quickly decay and only the velocity os- 
cillations remain. This explains why the power haloes are seen 
only in the velocity power maps. 

Thus, our model suggests that the high-frequency acous- 
tic oscillations should be trapped around the magnetic network 
cores and the active regions, where acoustic cavities can be 
formed under small-scale cylindrical magnetic canopies. The de- 
duced oscillation periods are within 1 .5 - 2 min for cavities with 
granular sizes, but additionally the cavities should oscillate with 
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the photospheric cut-ofF frequency. However, the picture will be 
different far away from the network cores and active regions, 
where the cylindrical bipolar magnetic fields are probably re- 
placed by almost a horizontal magnetic canopy which has been 
used as a model for a long time (Evans & Roberts, 1990). In 
this case, the higher frequency oscillations propagate upwards 
(only trapped cut off modes appear) and, consequently, power 
halos can not be formed. The oscillation at cut-off frequency 
should have a small amplitude there, as the rectangular verti- 
cal cavity does not have an oscillation mode with this frequency. 
The wake oscillating at cut-off frequency should be observed 
there, but with a smaller amplitude as almost the whole ini- 
tial energy is carried away by an initial pulse or a wave train. 
Therefore, it is the cylindrical structure of the field-free regions 
(see Fig. 1), which helps to trap acoustic waves as the oscilla- 
tions occur partly along the 0-direction. This explains why the 
"power halos" and "magnetic shadows" are observed only near 
the quiet-Sun chromospheric magnetic network cores (Mcintosh 
and Judge, 2001; Krijger et al., 2001; Vecchio et al., 2007). 

It must be noted that intensity halos are visible in the mag- 
netic canopy regions in some observations (Moretti et al. 2007). 
However, this is not inconsistent with our model. As mentioned 
before, our calculation does not include the exact solution in the 
overlying magnetic canopy region due to mathematical difficul- 
ties. On the other hand, the free oscillations of the cavity/canopy 
interface may excite oscillations in the overlying region as ob- 
served by Moretti et al. (2007). Additionally, the oscillations of 
the cavity may excite transverse waves in overlying magnetic 
canopy, which can be checked by observations of the magnetic 
field oscillation in the power halo regions. 
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Appendix A: Derivation of the velocity components 

Integration of Eqs. (37) and (38) gives: 



= e"'"-/^- su.4>^^ J- f ""'^':J"^^^ dt', (A.1) 



2Apo V n 




2Apo ynj 



f3/2 



where t' = (Oact. 

Integral in Eqs. (A.l), (A.2) can be solved as 



/cos(f' - 2 
?372 



= COS 



3n 



cos 



■ sm ■ 



37r r cos /' , 3;r r sin f , , 

2^[2nS[^2t'ln) + Fi 
2 V2^C ( yjlt'ln) + F2 



2 cos t' 



4 




37r 


2 sin t' 


T 





2 



37T . 37T 

COS — COS t +sm — sin f 
4 4 



yf2 

+2^/^[s{^|2FJ?) + c{^|2FJ?j] + —(F2-Fl) 
cos(f-3n/4) . 2^[S{^^C{^] 



= -2 - + 

V2 



jn/2 

■(F2-F1), 



(A.3) 



where Fi, F2 are the integration constants, and 5, C is the 
Fresnel sine and cosine integral functions, respectively. 

Fresnel functions tend to 1/2 for large arguments, therefore 
we have 

Ur, should tend to zero when f — > 00, therefore 

— (F2-Fi) = -2^^, 

which gives 

,cos(?' - 3;r/4) 



(A.5) 



/ = -2- 



f'l/2 



(A.6) 



Substitution of the expression (A.6) into Eqs. (A.l), (A.2) 
gives Eqs. (39) and (40). 
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